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ABSTRACT 

The orbital period of the binary pulsar PSR J0045-7319, which is located in our 
neighbouring galaxy the Small Magellanic Cloud (SMC), appears to be decreasing on a 
timescale of half a million year. This timescale is more than two orders of magnitude 
smaller than what is expected from the standard theory of tidal dissipation. Kumar and 
Quataert (1997a) proposed that this rapid evolution can be understood provided that 
the neutron star's companion, a main sequence B-star, has set up significant differential 
rotation. They showed that the spin synchronization time for the B-star is similar to 
the orbit circularization time, whereas the time to synchronize the surface rotation is 
much shorter, and thus significant differential rotation in the star is indeed expected. 
However, their calculation did not include the various processes that can redistribute 
angular momentum in the star, possibly forcing it into solid body rotation; in that case 
the dissipation of the tide would not be enhanced. The goal of this paper is to include 
the redistribution of angular momentum in the B-star due to meridional circulation and 
shear stresses and to calculate the resulting rotation profile as a function of time. We 
find that although angular momentum redistribution is important, the B-star continues 
to have sufficient differential rotation so that tidal waves are entirely absorbed as they 
arrive at the surface. The mechanism proposed by Kumar and Quataert to speed up the 
orbital evolution of the SMC binary pulsar should therefore work as suggested. 
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1. Introduction 

An interesting pulsar binary system was discovered 
in a systematic search for pulsars in the Small Magel- 
lanic Cloud (SMC) by McConncll ct al. (1991). This 
pulsar, PSR J0045-7319, has a spin period of 0.93 s 
and it is the only known pulsar in the SMC. Its com- 
panion is a main sequence B-star of mass ~ 8.8 Mq 
(Bell et al. 1995). The orbital period of this binary 
system is 51 days and its eccentricity is 0.808. Accord- 
ing to Kaspi et al. (1996), the orbital period of this 
system is decreasing on a time scale of ~ 5 x 10 5 years, 
whereas the standard theory for tidal dissipation pre- 
dicts an evolution time of order 10 9 years (much of 
our understanding of dynamical tide is due to a se- 
ries of papers by J. -P. Zahn, cf. Zahn 1975). The 
problem with the standard tidal theory, as applied to 
a rigidly rotating B-star of the SMC binary, is that 
only a small fraction of the tidal energy is dissipated 
in one orbital period. Kumar and Quataert (1997a, 
hereafter KQ) pointed out that if the entire tidal en- 
ergy were to be dissipated in about one orbital pe- 
riod then the resulting evolution of the orbit would 
be consistent with the observation. They suggested 
that this can be accomplished provided that the rota- 
tion in the interior of the B-star deviates significantly 
from the rigid body rotation. This is in fact what 
is predicted by the work of Goldreich and Nicholson 
(1989) who showed that the angular momentum car- 
ried by tidal waves is deposited first near the surface 
of the star where the waves are dissipated. Thus, the 
surface can be brought into synchronous rotation on 
a time scale short compared to the time it takes for 
the star as a whole to become synchronous. The time 
scale for the latter process for the SMC pulsar bi- 
nary system is shown to be comparable to the orbital 
circularization time (Lai 1996, Kumar and Quataert 
1997a) and thus the star is not expected to be rotat- 
ing synchronously; the observed apsidal motion of the 
star provides support for this picture (see Lai et al. 
1995, Kaspi et al. 1996). However, if angular mo- 
mentum is efficiently redistributed in stellar interiors, 
then a star can continue to rotate rigidly inspite of 
the tidal torque operating at the surface. The goal of 
this paper is to access the efficiency of angular mo- 
mentum redistribution and to calculate the expected 
magnitude of differential rotation for the B-star of the 
SMC binary system. In §2, we describe our models 
and the physical mechanisms considered in this study. 
The main results are summarized in §3. 



2. Angular velocity profile of a star due to 
tidal torque, meridional circulation and 
shear stresses 

We consider the evolution of the rotation profile 
of a main sequence B-star of 9 Mq. We consider the 
star to be rotating uniformly until such a time so that 
the radius of the B-star, result of the normal 

stellar evolution, becomes ~ 6 Rq. Fig. 1 shows the 
radius of the star and the size of its convective core 
as function of time. 

The energy in the dynamical tide depends on the 
rotation rate of star. The tidal energy in a B-star, for 
parameters corresponding to the SMC binary system, 
and angular rotation speed perpendicular to the or- 
bital plane of —1.6 x 10~ 5 rad s _1 (6 x 10~ 6 rad s~ 4 ), 
is about 6 x 10 40 erg (10 40 erg) (see KQ); the minus 
sign refers to retrograde rotation. The average angu- 
lar momentum luminosity associated with the tidal 
gravity wave for these two cases are 3 x 10 38 gm cm 2 
s~ 2 (10 38 gm cm 2 s~ 2 ) respectively. 

The dominant process for tidal wave absorption 
in early type stars is the radiative damping which 
arises when photons diffuse from regions that are com- 
pressed to regions that are expanding. Since the pho- 
ton mean free path increases rapidly with distance 
from the center of the star, most of the dissipation 
of the wave occurs close to the stellar surface. The 
fraction of wave momentum luminosity absorbed can 
be calculated using /damp = exp [T(r M )] where 




(1) 

I = 2 for the quadrupole tide, K = c/pn is the ther- 
mal diffusivity, n is the lower turning point defined 
by N(ri) — u>, r u is the upper turning point defined 
by r u — 6ic s /uj (n < r u ), c s is the sound speed, c 
is the speed of light, and A is the Brunt- Vaisala, fre- 
quency which can be written as a sum of two parts, 
one of which is depends to the thermal gradient (V) 
and the other on the gradient of the mean molecular 
weight (V M ) i.e. 

A 2 = A 2 + A 2 = -g-(V ad - V) + gv„ (2) 

(see Unno et al. 1989, or Zahn et al. 1997 for details). 
The local luminosity of the wave is then given by 

Cj(u, £, r) = Cj(u, £, n) exp[-r(w, £, r)]. (3) 
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If the star is rotating as a solid body, only the 
waves of frequency less than about 1.2 x 10~ 5 rad s _1 
are completely absorbed as they reach the surface. 
For higher frequencies, the fraction of the wave flux 
absorbed decreases as ~ a; -5 . However, differential 
rotation modifies that picture. Then, the frequency 
for the quadrupole tidal wave in the local rest frame 
of the fluid is uj(r) = u> t — 2Q(r), where f2(r) is the 
angular speed of the star, uj t « 2f2 p is the tidal fre- 
quency as seen in the inertial frame, and £l p is the 
orbital angular speed of the star at periastron. In a 
differentially rotating star, the local tidal wave fre- 
quency decreases outward and the wave dissipation 
rises rapidly. As waves deposit their angular momen- 
tum near the surface, the "synchronized front" starts 
at the surface and progresses inwards. However, one 
must also consider the role of the redistribution of 
angular momentum throughout the star as a result of 
meridional circulation and shear stress. These pro- 
cesses are calculated using the procedure described in 
Talon et al. (1997). 

The evolution of the angular momentum in the 
model is governed by 
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where p is the density, v v is the vertical (turbulent) 
viscosity and U is the amplitude of the vertical cir- 
culation speed. The latter is calculated considering 
differential rotation in radius, as described in Zahn 
(1992) and in Matias et al. (1997), neglecting the role 
of the mean molecular weight gradients. The gradient 
of mean molecular weight would somewhat increase 
the timescale for the redistribution of momentum by 
the meridional circulation in the central part of the 
star. However, it would influence only marginally the 
results for the outer portion, which is the main ob- 
ject of this paper. The turbulent viscosity takes into 
account the weakening effect of the thermal diffusiv- 
ity as first described by Townsend (1958). The details 
concerning that viscosity are not important here since 
the circulation, being an advective process, dominates 
the transport of momentum almost everywhere. 

The local time scale that determines the redistri- 
bution of angular momentum due to meridional cir- 
culation is the Eddington-Sweet time scale. Fig. 2 
shows this time scale as a function of r in the radia- 
tive part of a 9Mq main sequence star of age 19.7 



Myrs. Since the Eddington- Sweet timescale is much 
shorter in the outer portion of the star, we expect 
the outer region to be rotating almost as a solid body 
whereas differential rotation can persist in the inner 
partj^. 

Equation is solved to determine the rotation 
profile as a function of time. For the model with 
mild retrograde solid body rotation, the upper turn- 
ing point of the tidal wave is located at a radius of 
0.85i?* and only about 1.0 % of its angular momen- 
tum luminosity is deposited during one orbital period. 
In this case, meridional circulation is efficient in re- 
distributing the angular momentum in the interior of 
the star. In fact, as long as less than ~ 30 % of 
the angular momentum luminosity of the tidal wave 
is deposited in the star, meridional circulation car- 
ries away almost all of this angular momentum and 
spreads it over much of the star, and thus, little differ- 
ential rotation can be maintained. When wave damp- 
ing is inefficient, the result of the rotation evolution 
calculation is dependent on the uncertain initial ro- 
tation profile of the star, and the time dependence of 
the angular momentum luminosity of the tide. 

For these reasons we have carried out a restricted 
evolution calculation. We start with a state of dif- 
ferential rotation in the star so that a substantial 
fraction of the angular momentum luminosity of the 
tidal wave is deposited in the outer part of the star. 
We then calculate how the meridional circulation, and 
shear turbulence modify the rotation profile. In par- 
ticular we wish to find an initial rotation profile lead- 
ing to an increase of the differential rotation; a de- 
crease in the differential rotation would decrease the 
rate of energy and angular momentum deposit and 
thus increase the time scale for the orbital evolu- 
tion. To this end our calculation was started with 
the rotation profile shown in Fig. 3a (at 20 Myrs). 
In that case, the rate of angular momentum deposit 
was about 50 % of the luminosity. The initial rota- 
tion profile for prograde rotation was solid-body, as 
shown in Fig. 3b; in this case, about 80 % of the an- 
gular momentum luminosity of the wave is deposited 
near the upper turning point (at ~ 0.96 R*). The 
subsequent evolution of the rotation profile in each 
of these two cases is shown in Fig. 3; note that the 
difference between the rotation rate in the core and 



1 Actually, meridionnal circulation does not tend to produce solid 
body rotation but rather a state of mild differential rotation 
with O CO rc /^surface ~ 1.2 for a 9 Mq star (see Talon et al. 
1997 for more details). 
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the surface is increasing with time until the surface 
is nearly synchronized. In Fig. 4, wc show the wave 
damping as a function of r for every rotation pro- 
file shown in Fig. 3. Note that in all of these cases 
the tidal wave energy is completely dissipated as the 
wave approaches the stellar surface, and thus the an- 
gular momentum transported by the meridional circu- 
lation and shear induced turbulence does not modify 
the rate of change of the orbital energy for the SMC 
binary pulsar system PSR J0045-7319 estimated by 
Kumar and Quataert (1997). Thus, the time scale for 
the evolution of the orbital period calculated by these 
authors is unaffected by the redistribution of angular 
momentum in the star. 

We have also calculated the effect of the differential 
rotation induced turbulence on the tidal wave damp- 
ing. We find that it can be an important contributor 
to the wave damping at initial stages of evolution, 
when the surface of the star is rotating differentially, 
helping to enhance the fraction of tidal angular mo- 
mentum luminosity deposited in the star. 

We would like to point out one simplification made 
in our calculation. Even though the energy in the 
tidal waves is modified as the rotation profile of the 
star evolves we took it to be constant. However, 
the error made is not large since the tidal waves 
are mostly excited near the boundary of the convec- 
tive core and radiative envelope where the rotation is 
nearly constant. 

We note that the contribution to the wave damping 
from the region near the bottom of the radiative zone 
is negligible inspite of the fact that the mean molec- 
ular weight has a large gradient there leading to a 
sharp increase of the Brunt- Vaisala, frequency. The 
reason is simply that the photon mean free path here 
is also very short and it takes a long time for photons 
to diffuse across a wavelength. Even if the dissipa- 
tion in that region is about 5 times greater than just 
outside of it, the contribution to the total damping 
remains negligible. 

3. Conclusion 

We have calculated the evolution of the rotation 
profile of a 9 M Q main sequence star which is sub- 
ject to a tidal torque from its companion neutron star 
and in which angular momentum is redistributed by 
meridional circulation and shear stresses. We calcu- 
lated models corresponding to two different rotation 
states: one model has a retrograde initial rotation of 



— 1.6 x 1CP 5 rad s _1 , and the other has a prograde 
initial rotation of 6 x 10~ 6 rad s _1 . The magnitude 
of the tidal torque was taken to be 3 x 10 38 gm cm 2 
s~ 2 for the retrograde rotation and 10 38 gm cm 2 s~ 2 
for the prograde rotation. The location of the dissipa- 
tion of the momemtum luminosity was calculated as 
a function of the rotation profile. The meridional cir- 
culation and shear stresses were calculated following 
the method described in Talon et al. (1997). We find 
that the fractional difference between the surface and 
the interior rotation rate in the star is a factor of a 
few. With this differential rotation, the frequency of 
tidal waves in the rest frame of the star is low enough 
that they are completely absorbed before they reach 
the upper turning point. For retrograde rotation of 
frequency about 6 x 10 -6 rad s _1 or greater and a 
significant amount of differential rotation, the energy 
dissipated per orbit is about 10 41 erg, which leads to 
orbit evolution time of <~ 5 x 10 5 years, consistent 
with the observations. 
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Fig. 1. — The radius of the B-star as a function of 
time is shown in the top panel, and the lower panel 
shows the size of the convective core (fractional mass) 
as the star evolves. The rapid change in the stellar 
radius and the mass of the convective core at about 
24.5 Myrs is the because the star is evolving off the 
main sequence. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v4.0. 
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Fig. 3. — Rotation profile of the star at several times 
(in millions of years). The top panel shows the case 
where the initial rotation rate of the star was retro- 
grade with angular speed of 1.6 x 10~ 5 rad s^ 1 , and 
the lower panel corresponds to the progradc rotation 
speed of 6 x 1CP 6 rad s _1 . 
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Fig. 4. — Relative magnitude of the momentum lumi- 
nosity C(r)/ ' C(r c ) = cxp(— t) as a function of time for 
the two different cases of initial rotation of the B-star: 
retrograde rotation with angular speed 1.6 x 10~ 5 rad 
s _1 , and prograde rotation rate of 6 x 10~ 6 rad s _1 . 
Note that with increasing time, as the rotation profile 
evolves and the thinkness of the synchronously rotat- 
ing shell near the surface increases, the wave dissipa- 
tion occurs at smaller r. 
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